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We suggest a method of derivation of the long-wave action of the model spin 
Hamiltonians using the non-linear partial differential equations of motions of the 
individual spins. According to the Vainberg's theorem the set of these equations are 
(formal) potential if the symmetry analysis for the Frechet derivatives of the system 
is true. The case of Heisenberg (anti)ferromagnets is considered. It is shown the 
functional whose stationary points are described by the equations coincides with the 
long- wave action and includes the non-trivial topological term (Berry phase). 
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The most of rigorous results in the dynamics of a one- or two-dimensional isotropic 
Heisenberg antiferromagnet due to a presence of non-trivial topological term in the action 
(Berry phase) has been obtained in a coherent state path integral formulation 1 . This ap- 
proach starts with the derivation a path integral representation for the partition function of 
an individual spin Z = Tr exp(—/3H(S)) where the Hamiltonian H(S) is a linear function 
of the spin operator S. Further, this result is generalized for a lattice problem using the 
fact the Hilbert space of the many spin system is the direct product of the one-spin Hilbert 
spaces at each site. In path-integral formulation is convenient to use an overcomplete basis 
of spin coherent states \N) labeled by the points N on the surface of the unit sphere with 
the expectation value of the spin (N\S\N) = SN. The exponential in the expression for 
Z is written as the Trotter product of a large number of exponentials each evolving the 
system over an infinitisimal imaginary time interval At and using the resolution of unity 
1 = / y^-\N)(N\ between the intervals (the integral is over the unit sphere). This procedure 
leads to the following representation for the partition function 





" P 




d 

dr 


j DN(t) exp 


jir 


(n(t) 










N(t)\-H(3N{t)) 



(1) 



where H(SN(t)) is the Hamiltonian in spin coherent representation. The first term in the 
action is the Berry phase term Sb and represents the overlap between the coherent states at 
two infinitesimally separated times. The manipulation Sb into a physically more transparent 
form may be found, for example, in RefA The result is 



S 



B 



P 1 

iS [ dr [ du \ N 



o 



o 



dN ON 
du dr 



(2) 



where u is a dummy integration variable, N(t,u)) is defined by 



N(t,u)) = exp (u (z(t)S + - z*(t)S~)) \N = (0,0,1) 



with the properties N(t,u = 1)) = N(r)) , N(t,u = 0) 



(3) 

N = (0, 0, 1)) and the com- 
plex number z is given in spherical coordinates z = — ^exp(— i<f>). The Sb presents the 
oriented area contained within the closed loop defined by N(t), < r < (3, N(0) = N(j3) 
on the spin sphere. 

In this Letter, it will be demonstrated how the long-wave action for spin Hamiltonians can 
be derived by using the non-linear differential equations of motion of spins via the solving 
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of the inverse problem of the calculus of variations: that is, finding the functional whose 
stationary points are described by the descriptive equations. 

As is well known the most fundamental question is whether the functional exists for a 
given operator. In the Ref.-^ Atherton and Homsy have suggested a stringent and elegant 
formalism based on Vainberg's theorem^ to derive the operational formulas for the case of 
an arbitrary number of nonlinear differential equations in an arbitrary number of indepen- 
dent variables and of arbitrary order to determine whether given differential operators are 
potential. 

As the simplest illustration of the method we consider the case of a Heisenberg ferro- 
magnet with the spin Hamiltonian H = — \ Y JimSiS m where the Si are spin S quantum 

(im) 

spin operators on the ith sites of a n-dimensional hypercubic lattice. The exchange constant 
Ji m > is non-zero just for the nearest neighbors. The non-linear differential equations 
describing the dynamics of Heisenberg ferromagnet can be obtained by taking the diagonal 
matrix elements of the equation of motion for the raising operator %h—^~ = [S^, H] of the 

N 

ith spin in spin-coherent representation \Q) = T7 |$i> > where < 0j < it and < 0j < 2ty 
parametrize the spin states on the unit sphere^. In continuum approximation the coupled 
c-number equations for the real variables 6 and <fi depending on the time t £ [ii,^] an d 
space coordinates r £ 3? n can be written as following (hereafter the lattice constant is unit) 

= —S sin 9-^- H — — ( A© — cos © sin (V0) 
at ri \ v ' 

Q = S— + ^- (2 cos © (V0V©) +sin©A0). (4) 

These equations present a vector of differential operators BP(zk) = (j,k = 1,2 and 
zi = 0, z 2 = 4>) 

R 1 = -S sin ©^ + ^ (a& - cos sin (V0 X 2 



dt h 

f)f\ TC2 

R 2 = S sin 0— + — (2 sin cos (V0V©) + sin 2 0A0) . (5) 

The second equation in (JIJ is multiplied by the factor sin that is necessary for the further 
symmetry analysis. The inverse problem is solved if an operator can be proven to be a 
potential operator, i.e. an analysis of the Frechet derivatives is symmetrical (Vainberg's 
theorem). According to the procedure in^ we evaluate the "two-tensor" of Frechet derivatives 
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in the directions u k 

R\ = lim - (R?(z k + eu k ) - R j (z k )) , (6) 



e 

R\ = - Ul ScosQ^ + ^ ((A Wl ) - Wl cos 20 (V0) 2 ) , (7) 
R] 2 = -S sin 6-^ - — sin 20 (V0Vw 2 ) , (8) 

R\ = 5 sin0— - + M1.SCOS0 — 
- 1 dt dt 

JS 2 

+^ (M!sin2eA0 + sin26 (V«iV0) + 2«i cos 26 (V0V0)) , (9) 



R\ = If- S in 2 6Aw 2 + — - sin 26 (Vu 2 VO) . (10) 
n a v ' 

The symmetry requirement is satisfied if / ipjR J k u k dV = j UjR 3 k ip k dV for every ipj, u k in the 

*2 

range of the i?- 7 . The integral is taken over all space coordinates and time / dV = / di / dr. 
The respective ordering of the vectors R? and z k is very important. It determines the form 
of the variational integral and required symmetry conditions. The symmetry test yields to 
the following boundary conditions 

J faRfadV = J UiR^dV => ^(^iVui - u^ip^dS = 0, (11) 

s 

J ^ x R) 2 u 2 dV = J u 2 R\-i\) X dV =>• j> sin 2G^ 1 u 2 V(f)dS = 0, -01^ sin 6|£ = 0, (12) 

s 

J %p 2 R 2 2 u 2 dV = J u 2 R 2 2 %p 2 dV =>- y sin 2 6(^ 2 Vw 2 - u 2 V^ 2 )dS = 0. (13) 

s 

The variational principle for the differential system can be constructed via the formula 
i . i 

/ dV j d\ ZiR % (\zi) where / d\ represents integration over the scalar variable A. It results 

o o 
in the following 



F 1 = J dV J dX GR^XG) = JdV J dX(-S sin(A6)) 



o 

JS 2 



d(XG) d<p 
dX ~dt 



oh ■ J dV (sin 2 6 (V0) 2 + (V6) 2 ) + d(0), (14) 
F 2 = J dV J dX (j)R 2 (X(f)) = J dV iS sin 6—0 - — sin 2 6 (V</>) j+C 2 (6). (15) 
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The C12 are some functions which are determined by the condition Fi = F 2 = F(Q, <fi). The 
using of the standard parametrization 

5(A) = SN(X), N(X) = (sin(A6) cos(0), sin(A6) sin(0), cos(AG)) , (16) 

— * — * —* 

where N(l) = N(t) is the physical field and N(0) = (0, 0, 1) is the north pole of the physical 
sphere yields to the Berry phase term in 

' d{XQ) dcf) 



S B = JdV J dA(-Ssin(A0))- 

(17) 



d\ dt 



tl 1 
-S I df I dt I dX\N 



dN ON 
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h o \ L 

The comparison of Eqs. ()14[ IT3j) gives the condition (f){t\) = 0(t 2 ). In contrary to a path 
integral formulation the claim N(ti) = iV(t 2 ) arises as a necessary condition of potentiality 

of the differential operator. The following surface integrals must be zero 5*1=/ dSOVO 

s 

and S*2 = § dScf) sin 2 OV0. The unifying of Eqs. lfTH ITo]) using the divergence theorem gives 
the final form of the functional 

1 * 2 

F = S B --J dfj dt (Q\H\Q), (18) 

where the part 

(Q \H\ Q) = -JS 2 (n-- sin 2 6 (Vcpf - l - (V0) 2 ) (19) 

is the density of the energy of the system. The expression (fT%|) presents the long-vawe action 
(Euclidean) if to replace t — > it (t G [0, 1/T]) and multiply the functional by —%. 

The same analysis can be made for a generic Heisenberg two-sublattice antiferromagnet 
( J < 0) with the non-linear differential equations forming a vector of differential operators 

Brh 7S* 2 / 2 2\ 

= -S'sine 1 -^ + — -coseisine 2 cos(0i-0 2 ) (n- (V9 2 ) - (V0 2 ) j 

JS 2 JS 2 
+ — — cos Oi sin 2 sin(0i — 2 ) A 2 H — ;— cos 61 cos 2 cos(</>i — 2 ) A 6 2 
n, n 

US 2 -> -> JS 2 

H — - — cos 61 cos 9 2 sin(0! — 2 ) ( V0 2 V6 2 ) H — — sin G>i sin 9 2 A 9 2 
k v ' k 

— sine lC ose 2 (n - (V0 2 ) ) = 0, (20) 
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R 2 is the analogous equation with substitution of the lower indexes 1^2, 
09 1 JS 2 



S sin 9i- 



dt 



+ 



sin9isin9 2 sin(0 2 — <f>i) in — (V9 2 ) — (V0 2 ) 



JS 2 JS 2 
+ sin 9i cos 9 2 sin(0 2 — 0i) A 9 2 H — — sin 9i sin 9 2 cos(0i — 2 ) A 2 



+ 



h 

2JS 2 



h 



sin 9i cos 9 2 cos(0i — 2 ) (V9 2 V0 2 ) = 



(21) 



and the same equation i? 4 with 1 ^ 2. Here, one must suppose z\ = Qi, z 2 = Q 2 ,z 3 = 
01,^4 = 2 where the lower indexes in 9 and denote the sublattices. The symmetry 
test / UiR\u 3 dV = J u 3 R\uidVand J u 2 R 2 A u 4 dV = J u 4 R 4 2 u 2 dV gives the time boundary 

conditions sin9i(ti) = sin9i(i 2 ) and sin9 2 (ti) = sin9 2 (t 2 ), correspondingly. The rest 

symmetry conditions determine the surface integrals which must be zero if the differential 

operator is potential. The construction of the variational principle gives 

i i 
Fi = J dV J dA9 lJ R 1 (A9 1 ) = -S J dV J dX | ,Y 
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J rfVsin9i (sin9 2 cos(0i - 2 ) (n - (V0 2 )~ - (V9 2 

+ cos 9 2 cos(0i — 2 ) A9 2 + sin 9 2 sin(0i — <f> 2 ) A0 2 
+2 cos 9 2 sin(0i - 2 ) (V9 2 • V0 2 )) 

dVcosQ 1 (cos9 2 (n-(V9 2 ) 2 j -sin9 2 A9 2 ) + d^, 9 2 , 2 ) 



(22) 



and 



l l 
dV J dA0 lJ R 3 (A0 1 ) = -S J dV J d\\ X 



dm dm 

x 



d\ 

2 



dt 



J S f { { 2 2 

+ — J dV sin 9i ( sin 9 2 cos(0i - 2 ) (n - (V0 2 ) - (V9 2 ) 

+ cos 9 2 cos(0i — 4> 2 ) A9 2 + sin 9 2 sin(0! — 2 ) A0 2 
+2 cos 9 2 sin(0i - 2 ) (V9 2 • V0 2 )) + C 2 (Oi, 9 2 , 2 ), 



(23) 



where C± i2 are some functions. The comparison of both forms yields to F(9i,0i) + 
Cs(9 2 , 2 ) = Fi = F 2 if the following conditions are fulfiled 



Ci(0i, 9 2 , 02 ) = C , 3 (9 2 , 2 ), 



(24) 
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c 2 (e 1 ,e 2 ,0 2 ) = c 3 (e 2 ,0 2 ) 



+ 



JS 2 



J dVcos& 1 (ros© 2 (n- (V6 2 )^ -sine 2 A0 2 



(25) 



The similar expression may be obtained for -F(0 2 ,0 2 ). It equals to F(Oi,0i) if to replace 
1^2. The final form for the variational principle may be obtained via an integration by 
parts using the divergence theorem and the boundary conditions found at the symmetry 
analysis. The result is 



F = -SY, df dt dXlNi 



dNi dNi 



1-2 



^ / df I dt(tt\H\tt) . (26) 



As in the previous case the expression coincides with the action and includes the Berry 
phase terms from each sublattices. The density of the energy is given as following 

(Vl\H\Vl) = JS' 2 (— n sin 0! sin 2 cos(0i — 2 ) — n cos 0i cos 2 
+ COS ©i cos 2 cos(0i — 2 ) (V©i V0 2 ; 
— sin ©i cos 2 sin(0i — 2 ) (V0i V0 2 ^ 
+ cos ©i sin 2 sin(0i — <f> 2 ) (V©i V0 2 ^ 



+ sin ©i sin @ 2 cos(0i — (j. 
+ sin ©i sin 2 (v©iV©; 



V0lV0 2 



(27) 



In conclusion, the method of derivation of the long-wave action for spin Hamiltonians is sug- 
gested via the solving of the inverse problem of calculus variations. The case of a Heisenberg 
(anti)ferromagnet is considered. 
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